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ABSTRACT 


Recent measurements of the long-term frequency stability 
of five commercial cesium beam frequency standards relative 
to the NASA NX-1 experimental hydrogen maser have yielded 
a fractional frequency difference between the two of 
(6.77 ± 1.06)X 10“^. Since the frequency of the NX-1 
maser was synthesized to the value for the hydrogen maser 
frequency given by Vessot et al. in 1966, this yields a new 

value for the hydrogen maser frequency of 1,420,405,751.7767 ± .0024 
Hertz. This result is within the stated error limits of several 
other recent measurements of the hydrogen maser frequency. 
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FREQUENCY COMPARISON OF 
FIVE COMMERCIAL CESIUM STANDARDS 
WITH A NASA EXPERIMENTAL HYDROGEN MASER 


INTRODUCTION 

Recent measurements of the long-term frequency stability of five commercial 
cesium beam frequency standards* relative to a NASA experimental hydrogen 
maser, 2 constructed at GSFC and designated NX-1, have made possible a new 
comparison of these two types of atomic frequency standards. The result of 
this comparison is a value for the hydrogen maser frequency that is 6. 77 
parts in 10*2 lower than that given by Vessot et al. 2 in 1966. 

The frequency comparison of the atomic standards is based on continuous 
phase difference measurements, sampled hourly, over the 123 days between 
May 7 and September 7, 1968. The 2952 hourly phase difference data thus 
obtained were punched on cards for computer analysis. Frequency data 
were obtained by taking the time derivative of successive five-data-point 
least-square-fits of the phase data, using Koenig’s approximation formula. 4 
The frequency data were then reduced to yield the mean fractional 
frequency difference between each cesium standard and the hydrogen maser. 

These figures were finally corrected so as to obtain the results in terms 
of the A-l second. 

DATA REDUCTION 

To obtain the results shown in Table 1 , the frequency data were least- 
square fitted (LSF) to a linear equation for frequency as a function of time. 

The fact that the total frequency change over the entire period of measurement 
due to the time dependent term was less than 1 * indicates that the cesium beam 
standards exhibit no systematic linear frequency drift. However, in calculating 
the mean frequency difference, it was convenient to take the sum of the LSF fre- 
quency intercept at t=0 and one-half of the total frequency change due to the LSF 
time dependent term (these two components are given in Table 1). That this 
procedure correctly yields the mean of the frequency difference data can readily 
be proved. (See Appendix ) . 

CORRECTIONS 

The frequency output of the NX-1 hydrogen maser is synthesized to the 
value given by Vessot et al. 2 in 1966, corrected for differences in bulb shape, 
second order Doppler shift, and applied magnetic field. To achieve frequency 
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Table 1 


Calculated Mean Fractional Frequency Difference Between 
Cs Standards and Hydrogen Maser, from Phase Measurements (in parts per 10 ^) 


Cs Units 

(a) 

(U„ 

(b) 

•W).. 

T 

2 

J£l ! 

H) ! 
1 

<d) 

a 

1 

6.63 

-.046 


l 

6.58 | 

1.07 

2 

6.32 

-.281 


1 

6.04 1 

1 

1.48 

3 

5.97 

-.047 


5.92 | 

0.90 

4 

3.12 

-.758 


1 

2.36 1 

1.21 

5 

4.25 

-.556 


1 

3.69 l 
L 

1.49 


° rms = 1.25*10“^ (RMS for all five Cs units) 


(a) 



is the intercept of an LSF linear equation for frequency as a 
function of time 



is the average frequency change during the total period of 
measurement, T = 2952 hours. 

is the mean fractional frequency difference between each 
cesium standard and the NASA hydrogen maser, NX-1. 


(d) ® 


is the standard deviation of — from its mean. 


A_f 
f * 


outputs which are in terms of universal time, this synthesized NX-1 output and 
the output of each of the cesium standards are offset by an amount which is 
nominally -300xi0"10, the UT-2 offset for 1968, relative to the defined cesium 
resonance frequency based on the A-l second. 5 

In the NX-1 this offset is maintained at exactly -300 x io -10 . In the cesium 
standards, however, the amount of the offset was found to be slightly larger. 
Hence the measured fractional frequency difference between the cesium 
standards and the hydrogen maser must be corrected for the difference between 
the actual offset in the cesium standards and -300 x 10" 1®. 
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The offset in the cesium standards is maintained by an internal electronic 
synthesizer and a magnetic field fine tuner. The first two columns of Table 2 
show the contribution to the offset of each of these sources. The third 
column shows the effect on the offset of the magnetic field drift which was 
observed during the measurements. (The approximately linear behavior of 
this drift with time, shown in Figure 1, was used to compute its contribution 
to the offset. ) The remaining columns of Table 2 show the total offset of each 
cesium standard and the appropriate correction to account for the difference 
between this offset and -300 * 10"* 0 . 


Table 2 

Magnetic Field and Synthesizer Offsets in 
Cesium Standards and Consequent Corrections. 


<The Off-Set Frequency of UT -2 for 1968 Relative to the A- 1 Time 
Scale is -300 * 10 "^ Exactly.) 


Cs Units 

(a) 

Synthesizer 
Offset 
(parts per 
10 10 ) 

(b) 

H Field 
Offset 
(parts 
per 
10 10 ) 

(c) 

<5H Field 
Offset 
(parts 
per 
10 12 ) 

Total Offset 
(parts per 
10 10 ) 

Correction to be 
applied to Fractional 
Frequency Difference 
(parts per 10 * 2 ) 

1 

-301.8722 

1.860 

0.544 

-300.0068 

0.68 

2 

-301.8722 

1.850 

0.419 

-300.0180 

1.80 

3 

-301.8722 

1.851 

0.251 

-300.0187 

1.87 

4 

-301.8722 

1.847 

-0.034 

-300.0260 

2.60 

5 

-301.8722 

1.849 

— 

-300.0230 

2.30 


(a) The synthesizer offset is furnished by the manufacturer. 

(b) The magnetic field offset is calculated from the Breit-Rabi formula, 

i. e. — = 427H 2 where H = -tt — 3 , i' being the Zeeman frequency. 

I 00U x lu 

(c) The applied magnetic field in the transition region between the Ramsey 
cavities was found to exhibit a small upward drift, as shown in Figure 1. 


3 












APPLIED MAGNETIC FIELD 


44.4 


44.2 


Af = -20 Hz 


44.4 



44.4 - 
44.2- 



Af = 130 Hz 


44.0 


10 20 30 

MAY 


_L 

10 


_J L 

20 30 

JUNE 


10 20 30 

JULY 


10 20 30 10 20 

AUG SEP 


Figure 1 — Magnetic field measurements of cesium standards. (Af is the total 
frequency change due to the magnetic fie Id variation over the 130 days of field 
measurement; the serial number (SN) and current control potentiometer dial 
reading (DR) are given for eoch unit.) 


Cs4 Vs HM 
(SN 139; DR 732) 

Cs3 Vs HM 
(SN 137; DR 736) 


Cs2 Vs HM 
(SN 138; DR 744) 

Cs t Vs HM 
(SN 136; DR 729) 



After applying the corrections discussed above, the total fractional 
frequency difference between each cesium standard and the hydrogen maser, 


/ A f \ 

V f / Cs-HM’ * S aS ®' ven * n Table 3. The mean fractio nal fr equency difference 

for the five units is 6.77 * 10" 12 . The deviation, <5 , of | — j j 


Cs-HM 


for each 


cesium standard from this mean is also given in Table 3, yielding a S of 

±1.11 x 10" 12 . 


Table 3 


Corrected Mean Fractional Frequency Difference between 
Cesium Standards and Hydrogen Maser, in terms of the 
Hydrogen Maser Frequency given by Vessot et al. (in parts per 10 12 ) 


Cs Units 

Af 

f 

M, 

\ f /Cs-HM 

6 

1 

6. 58 

0.68 

7.26 

0.49 

<*. 

6.04 

1.80 

7.84 

1.07 

3 

5.92 

1.87 

7.79 

1.02 

4 

2.36 

2.60 

4.96 

-1.81 

5 

3.69 

2.30 

5.99 

-0.78 

( 1 ) 

v Cs-HM 

8 RMS = 1,11 

x 10" 12 (Mean for all five Cs Units) 
x 10“ 12 (RMS for all five Cs Units) 



ESTIMATED ERRORS 

If the total error in the mean fractional frequency difference is taken to be 
the root-sum square of »rmS, from Table 1, and 6 rms. from Table 3, a 
total error of 1.66 * 10" * 2 is obtained. This is consistent with the result 
obtained by combining the estimated absolute errors of the Hewlett-Packard 
cesium standards (1.33 * 10" 12 ) and the NX-1 hydrogen maser (1.0 * 10" 12 ), 
which yields a total of 1.64 * 10"* 2 . These individual errors are compared 
in Table 4 with those of other similar units. 
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Table 4 


Errors in Cesium and Hydrogen Maser Standards 


7„pe of Standard 

(a) 

<r 

10-12 

Source 

Comment 

Cs (HP5060 A- 506 1 A) 

1.33 

H/P 6 

Based on factory test of 150 units relative 
to the company's cesium house 
standard, 3 «r = 4 y 10 ”* 

Cs (NBS III) 

1.1, 1.9 

NBS 3 * 7 

Under normal conditions, 3 <r = 

5 . 6 x 10“12; during the 1966 measure- 
ments with the Varian H-10 maser the 
NBS III received special care and its 
error was estimated at 1 = 1. 1 * 10 “ 12 

Cs (NRC III) 

.75 

NRC 8 

Based on reevaluation of NRC Cs III, 
2 <r - 1. 5 x 10-12 

H Maser (NASA NX- 1 ) 

1.0 

NASA 9 

Based on repeated comparisons of NASA 
NX-1 and Varian H-10 masers 

H Maser (NRC) 

0.9 

NRC 10 

Based on inter comparison of two 
hydrogen masers, HM-1 and HM -2 

H Maser (H-10) 

0.47 

H/P 3 

Based on Varian H-10 #4 


(a) It should be noted that the 1 cr errors given in this table are not directly comparable, since the Hewlett- 
Packard figure represents the dispersion of a large number of production cesium standards, while 
the others are estimates, based on measurements of physical parameters which are known to 
affect the frequency, with the errors computed in a different way in each case. 




CONCLUSIONS 


Our measured value for the fractional frequency difference between the 
commercial cesium standards and the NASA experimental maser is thus 
(6.77 ± 1.66) x 10“ 12 . Since the frequency of the NASA maser is synthesized 
to the 1966 measurement by Vessot et al. , our value for the hydrogen maser 
frequency is therefore fjjM = 1,420,405,751.7767 ± .0024 Hertz. 

If this value is compared with two other measurements also made in 1968, 
those of A. G. Mungall et al. 10 and C. Menoud et al. 11 as shown in Table 5, it 
is evident that these three independent measurements agree within their 
stated errors. 

It is also interesting to note that the three results shown in Table 5 fall within 
the error limits of the earlier results of Winkler 12 and of Peters et al. 13 . The 
discrepancy between these five measurements and the measurement of Vessot et 
al. indicates that further comparisons are desirable. 


Table 5 


Comparison of Recent Results of Hydrogen Maser Frequency 
Measurements (since 1968) with the Value Given 
by Vessot et al. (in 1966) 


(Vessot et al. value = 1,420,405,751.7864 ± .0017 Hertz) 


Source 

Hydrogen Maser Frequency 
1,420,405,751 + 

Identifications 

GSFC 

.7767 ± .0024 

H/P5060&5061 vs NASA NX-1 

NRC 

.7763 ±.0020 

NRC HM-2 vs NRC CSIIl( a) 

LSRH 

.7782 ± .0036 

LSRH H-2 vs Ehanches CS-03^ 


(a) 

NRC 

(b) 

LSRH 


National Research Council, Ottawa Canada 

Swiss Laboratory for Horological Research, Neuchatel, 
Switzerland 
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APPENDIX 

PROOF OF FORMULA FOR MEAN FRACTIONAL FREQUENCY DIFFERENCE 


As mentioned in the text, the frequency data obtained were least square 
fitted to an equation linear with time, i. e. 


Af . . 

T “ a o V- 


(Al) 


The mean fractional frequency difference, ( -^r) , was then computed using the 


formula 




(A2) 


where T = 2952 hours, the total period of measurement. (The terms a Q and 
j a^T are given in columns 1 and 2 of Table 1 of the text; the resulting 

is given in column 3. ) The purpose of this appendix is to prove that equation 

(A2) is correct for the mean of . 

The linear least square fit of |(xj, yj)| is given by 

' y = a Q + a x x (A3) 

where 


a o 


ly l 

yx. 

i 

lx.y. 

2 

2Xi 

n 

SX i 


2 

2X i 

PCi, 


Zx.sy. -sx.sx.y. 

2 2 

n £x. - (2X.) 


(A4) 
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and 


n 


sx i 

n 

Sx. 

1 

sx i 

2 

SX i 


nsx.y. -ix.iy. 

2 2 
niXj - (£x.) . 


Theorem. Let y = a Q + a^x be a linear LSF of |(x., y^)|.” ^ Then 


(A 5 ) 


- I y = a„ + a 
n i=i 0 1 


(ihU 


(A 6 ) 


Proof: 


a \ swiw! n s x i y, - 

*0 +a l(iT ZX i) = 5~7 \2 + 2 7 \2 

' ' »S*, -(I*,) "*V P 



P xy i ■ - x i P y i +PP y i' n 1 

Pi) 

f IVi 

^ 2 1 
n2x. - 1 

(p) 

I 2 


m — 


■*?-(*«,)* 

(n Iy i, 


n 2 x 


f - W 


i " 

- 2 y. 

n i-i 1 


(A 7 ) 
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Now consider the case of equally spaced abscissa values, i. e. where 


x i = Xf + (i - l)h, (A8) 

for 1< i <n. This permits the following simplification of the average of the 
x.'s: 

1 n i n 

- Z x = - £ [x + (i - l)h] 

n i n * 


1 f 

= - nx, 
n L 1 



, (n - l)h 

X l + “T— 


2 i x i + t x ! + (n " 1)h H 


x, + x 
1 n 


Substituting this result into (A6) we have 


1_ 

n 


n 

I 

i=l 




(A9) 


i. e. , for equally spaced abscissa values, the value of the LSF at the midpoint 

of the interval [ x, , x ] is the average of the y., 1 < i < n. 

In i 


In our case x = t and y = ~; hence, taking = 0 and x n = T and using the 
definition 



(A10) 


equation (A9) reduces to equation (A2), Q. E.D. 
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